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This  work  analyzes  the  CO-PrOx  reactor  design  as  a  component  of  the  CO  clean-up  system  of  the  ethanol 
processor  for  H2  production  applied  to  PEM  fuel  cells.  The  operating  conditions  of  the  processor  require 
compact  and  lightweight  pieces  of  equipment  and  efficient  operation  at  different  conditions.  An  egg¬ 
shell  catalyst  type  of  Pt/Al203  is  considered.  One-dimensional  heterogeneous  catalytic  reactor  model 
accounting  for  interfacial  gradients  is  used  to  optimize  the  PrOx  reactor.  Different  reactor  components 
are  added  gradually  to  illustrate  how  the  system  dimensions  and  configuration  change  after  optimiza¬ 
tion.  The  optimization  problem  determines  the  optimal  reactor  length,  reactor  diameter,  catalyst  particle 
diameter,  inlet  reactants  temperature  and  insulating  material  thickness  that  minimize  the  total  system 
volume.  On  these  model-based  results,  the  final  reactor  design  is  mainly  governed  by  the  presence  of 
hemispherical  heads  (distributor  and  collector).  Different  inlet  CO  compositions  and  power  generation 
targets  are  analyzed.  According  to  the  inlet  CO  level,  more  than  one  catalytic  stage  is  required  to  meet 
design  goals  and  fulfill  process  constraints.  The  model-based  reactor  optimization  of  the  pseudo-adiabatic 
operation  allows  obtaining  both  designs  for  reducing  volumes  and  optimal  operating  conditions  that 
allows  conventional  reactor  technology.  Afterwards,  simulation  runs  based  on  a  rigorous  one-dimensional 
heterogeneous  catalytic  reactor  model  accounting  for  intra-particle  gradients  are  performed  using  data 
obtained  from  model-based  optimization  results.  The  aim  of  these  simulations  is  to  verify  feasibility  of  the 
optimal  design  obtained  from  the  proposed  one-dimensional  heterogeneous  catalytic  reactor  model  with¬ 
out  intra-particle  gradients,  which  is  intended  to  approximate  an  egg-shell  catalyst  behavior.  The  present 
work  reflects  clearly  the  advantages  of  applying  mathematical  programming  techniques  to  optimize  both 
design  and  operation  conditions  of  the  PrOx  reactor. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

During  the  last  decade,  there  have  been  important  advances 
in  fuel  cells  technology.  Fuel  cells  are  being  developed  for  appli¬ 
cations  to  electrical  energy  generation  and  co-generation  systems 
(coupled  heat  and  power)  in  both  stationary  and  mobile  systems. 
Since  fuel  cells  are  high-efficiency  energy  converter  devices  and 
because  of  their  low  polluting  emission  levels,  they  become  more 
and  more  attractive  as  a  power  generation  alternative,  specially  in 
the  transportation  industry  [1,2]. 

Proton  exchange  membrane  fuel  cell  (PEMFC)  demands  a  free- 
CO  hydrogen  stream  for  operation.  Since  CO  is  adsorbed  on  the 
catalyst  surface  causing  catalyst  poisoning,  it  is  necessary  to  reach 
CO  concentrations  of  less  than  10  ppm  for  preventing  irreversible 
damage  and  to  facilitate  the  electrochemical  reaction  on  the  Pt  elec¬ 
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trode.  The  conditioning  of  the  gas  stream  generated  from  the  steam 
reforming  of  ethanol  is  partially  performed  by  the  water-gas-shift 
reaction  (WGS).  The  WGS  stage  in  the  reforming  systems  allows 
achieving  an  outlet  gas  composition  with  a  typical  CO  composition 
of  0.5-2%  (v/v).  Afterwards,  a  final  CO  reduction  is  needed. 

Physical  methods  (pressure  swing  adsorption  (PSA))  and  cat¬ 
alytic  methods  (methanation  and  preferential  oxidation)  have  been 
studied  for  selective  CO  removal  from  a  hydrogen-rich  stream. 
Among  these  different  hydrogen  purification  systems  the  CO  pref¬ 
erential  oxidation  is  preferred  for  small-scale  reforming  systems 
because  of  its  relative  simple  implementation,  lower  operating 
costs  and  minimal  hydrogen  loss  [3,4]. 

For  fuel  cell  applications,  a  compact,  efficient  and  reliable  fuel 
processor  is  desirable.  Process  synthesis  and  design  tasks  are 
similar  to  other  industrial  reforming  processes.  Flowever,  the  pro¬ 
duction  capacity  level  required  by  fuel  processors  for  vehicles  or 
other  similar  devices  is  lower  when  compared  to  industrial  pro¬ 
cesses.  Thus,  process  units  of  small  size  and  specific  designs  are 
required. 
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Nomenclature 

av  external  catalyst  surface  area  per  unitary  reactor 

volume  (cm-1) 

Ar  cross-reactor  area  (cm2) 

Q  component  concentration  along  the  reactor 

(mol  cm-3) 

Cpf  heat  capacity  (J  g”1  K-1 ) 

CS  l  component  concentration  inside  the  catalyst  pore 
(mol  cm-3) 

C| .  component  concentration  at  the  external  catalyst 
surface  (mol cm-3) 

Di  k  binary  diffusion  (cm2  s-1 ) 

Dp  catalyst  particle  diameter  (cm) 

Dt  reactor  tube  diameter  (cm) 

Difi  m  diffusion  coefficient  of  the  component  i  in  the  gas 
mixture  (cm2  s-1 ) 

eins  insulating  material  thickness  (cm) 

et  reactor  tube  thickness  (cm) 

Fi  molar  flow  of  component  i  (mol  s-1 ) 

G  mass  flow  per  unitary  reactor  area  (g  cm-2  s-1 ) 

ftf  energy  transfer  coefficient  in  the  film 

(Jem-2  s-1I<-1) 

A Hj  heat  of  reaction  (J  mol-1 ) 

kgj  mass  transfer  coefficient  in  the  film  (cm  s-1 ) 

Lt  reactor  length  (cm) 

Mi  molecular  weight  of  component  i  (g  mol-1 ) 

P  pressure  (atm) 

Qex  heat  exchanged  between  the  reactor  and  its  sur¬ 
roundings 

vj  reaction  rate  (mol  g-cat-1  s-1 ) 

If  fluid  temperature  (K) 

Tins  surface  insulating  material  temperature  (K) 

Ts  temperature  inside  the  catalyst  pore  (I<) 

Tf  temperature  at  the  catalyst  surface  (K) 

Tci  critical  temperature  of  the  i  component  (I<) 

Vtotal  volume  (cm3) 

VC/  critical  volume  of  the  component  i  (cm3  mol-1 ) 

Greek  letters 

ol\  stoichiometric  coefficient  of  component  i  in  reaction 

j 

sb  bed  porosity 

kf  gas  thermal  conductivity  (J  cm-1  s-1  K-1 ) 

/Xf  gas  viscosity  (g  cm-1  s-1 ) 
ph  bulk  catalyst  density  (g-cat  cm-3 ) 
pi  gas  phase  density  (g  cm-3 ) 
pp  catalyst  density  (g-cat  cm-3 ) 


The  aim  of  this  work  is  focused  on  investigating  the  PrOx 
reactor  as  a  process  component  of  a  fuel  processor  for  appli¬ 
cations  in  PEM-type  fuel  cells,  and  showing  how  the  reactor 
design  using  mathematical  programming  techniques  allows  com¬ 
puting  both  reduced  volumes  and  optimal  operation  conditions. 
Different  reactor  configurations,  reactor  size  and  relative  sizes 
of  the  reactor  components  (insulating  material,  reactor  tube  and 
reactive  bed)  can  be  evaluated;  moreover,  process  performance 
bottlenecks  and  opportunities  for  optimization  can  be  identified 
and  analyzed.  This  knowledge  can  be  used  for  process  improve¬ 
ments  such  as  lower  unit  weight  and  costs,  and  higher  global 
efficiency. 


2.  Preferential  oxidation  of  carbon  monoxide 

Carbon  monoxide  oxidation  to  obtain  carbon  dioxide  in  a 
hydrogen-rich  stream  was  first  investigated  in  the  ammonia  indus¬ 
try  in  the  60s.  As  a  result  of  those  former  works,  a  highly  active 
and  selective  platinum-based  catalyst  -  named  Selectoxo  -  was 
developed  [3].  The  main  reactions  involved  are:  carbon  monoxide 
oxidation  (1 )  and  hydrogen  oxidation  (2): 

CO  +  i02  C02  AH°98  =  -283.2  kjmol-1  desirable  (1) 

H2  +  202  -*■  H20  AH°gg  =  -241  kj  mol-1  undesirable  (2) 

Some  researchers  pointed  out  that,  depending  on  the  operation 
conditions  and  the  catalyst  used,  the  inverse  water-gas-shift  reac¬ 
tion  rWGS  (3)  and  the  methanation  reaction  (4)  can  also  occur,  but 
to  a  lower  extent  [4,5]. 

CO  +  H20  4*  C02  +  H2  AH%g  =  -42.2 kjmol"1  (3) 

CO  +  3H2  CH4  +  H20  AH°98  = -205.3  kj  mol-1  (4) 

Excessive  hydrogen  consumption  is  undesirable  since  the  sys¬ 
tem  efficiency  diminishes,  and  makes  it  difficult  to  control  the 
system  temperature.  Thus,  the  catalyst  has  to  facilitate  or  favor 
reaction  (1)  over  reaction  (2)  [6].  The  catalyst  would  ideally  have 
to  selectively  oxidize  about  10,000  ppm  (1%)  CO  to  less  than  5  ppm 
without  oxidizing  any  of  the  present  30-70%  H2.  For  example,  for 
a  molar  ratio  02/C0  =  2  and  using  a  Pt/alumina-based  catalyst,  a 
100  ppm  CO  concentration  is  obtained  at  the  reactor  outlet.  How¬ 
ever,  using  the  same  catalyst  and  for  a  02/CO  ratio  equal  to  3,  the 
CO  removal  is  complete  in  a  mixture  containing  1%  CO.  But  the  02 
excess  also  consumes  H2,  consequently  decreasing  the  processor 
efficiency.  In  addition,  if  a  large  excess  of  air  is  supplied,  the  tem¬ 
perature  control  becomes  difficult  due  to  the  exothermic  nature 
of  the  reaction  [7].  Then,  it  is  necessary  to  have  a  highly  selective 
catalyst  that  allows  operating  at  02  levels  close  to  the  stoichiomet¬ 
ric  one.  The  selection  of  a  highly  selective  catalyst,  reactor  design 
and  determination  of  operation  conditions  are  critical  aspects  that 
have  to  be  addressed.  Once  the  catalyst  to  be  used  is  defined,  tem¬ 
perature,  02/CO  ratio,  reactor  operation  mode  and  the  number  of 
stages  (i.e.  reactors  in-series)  constitute  key  decisions,  leading  to 
opportunities  for  process  optimization. 

The  efforts  for  improving  the  PrOx  reactor  performance  have 
been  focused  in  a  wide  spectrum  of  subjects,  including  the  devel¬ 
opment  of  more  active  catalysts  [3,5,6,8-15],  experimental  and 
theoretical  studies  to  formulate  kinetic  expressions  for  novel  cat¬ 
alysts  [8,12,15-18]  and  theoretical  studies  related  to  the  effect  of 
heat  transfer  resistance  on  PrOx  reaction  behavior  [16].  Research 
is  also  concerned  with  applying  and  evaluating  non-conventional 
reactors  such  as  microreactors  [17]  and  monolith  reactors  [18]. 

Several  catalysts  have  been  formulated  to  perform  the  carbon 
monoxide  preferential  oxidation,  most  of  them  based  on  noble  met¬ 
als.  Researchers  have  reported  the  use  of  supported  catalysts  based 
on  iridium  (Ir)  [3],  palladium  (Pd)  [3],  platinum  (Pt)  [3, 5, 6, 8, 9],  gold 
(Au)  [8,10],  rhodium  (Rh)  [3],  rutenium  (Ru)  [3,11,12]  and  platinum- 
based  alloys  [3]. 

In  order  to  reduce  costs  and  improve  selectivity,  recent  works 
have  been  focused  on  using  non-noble  metals.  Sedmak  et  al. 
[13]  developed  a  nanostructured  non-stoichiometric  Cu0.iCe0.gO2- 
Y  catalyst.  Zhou  et  al.  [14]  investigated  the  catalytic  properties  of  Ni, 
Co  and  Co-Ni  supported  on  activated  carbon.  Marino  et  al.  reported 
that  the  catalytic  system  CuO-Ce02  is  very  active  and  selective  for 
CO  oxidation  in  the  presence  of  a  large  amount  of  H2  [15].  These 
metals  have  demonstrated  to  have  a  high  activity  for  CO  oxidation 
in  hydrogen-rich  streams. 
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Platinum-based  catalysts,  supported  on  alumina  or  silica  are  the 
most  used  and  investigated  as  conventional  catalysts  for  the  CO 
preferential  oxidation.  Although  the  oxidation  rate  of  H2  is  sev¬ 
eral  magnitude  orders  higher  than  CO  on  a  Pt  surface,  and  that  its 
concentration  in  the  gaseous  stream  entering  the  reactor  is  two 
magnitude  orders  greater  than  CO,  the  preferential  oxidation  of 
CO  is  possible  due  to  the  relative  values  of  the  adsorption  heats 
of  the  substances  involved  on  a  Pt  surface  [19].  Therefore,  all  these 
catalytic  systems  exhibit  a  range  or  “window”  of  operation  tem¬ 
perature.  Whereas  on  the  one  hand  sufficiently  high  operation 
temperatures  are  desirable  to  increase  the  reaction  rate  (1)  and 
to  diminish  consequently  the  reactor  size,  on  the  other  hand  it  is 
convenient  to  operate  at  sufficiently  low  temperatures  to  avoid  H2 
oxidation.  The  temperature  range  in  which  activity  and  selectiv¬ 
ity  requirements  are  satisfied  is  so  narrow  that  the  catalyst  cannot 
admit  small  system  disturbances.  Then,  the  previously  mentioned 
catalyst  types  require  a  precise  and  complex  temperature  control 
system  [20]. 

Zalc  and  Loffler  [21]  also  pointed  out  that  reactor  temperature 
control  is  crucial  due  to  the  reaction  kinetics  and  the  narrow  tem¬ 
perature  operation  range  to  achieve  high  selectivity  for  the  desired 
reaction  (1 ).  A  single-stage  unit  operating  in  adiabatic  mode  would 
yield  poor  system  performance.  Thus,  an  optimal  reactor  design  has 
to  include  an  effective  temperature  control.  To  accomplish  this,  they 
highlight  that  the  following  non-excluding  alternatives  can  be  con¬ 
sidered:  supporting  catalyst  on  a  heat  exchanger,  or  using  air  staged 
injection  or  water  injection.  They  also  estimated,  based  on  simula¬ 
tions,  a  catalytic  bed  volume  of  approximately  3.5  L  for  a  hydrogen 
production  of  around  31  mol  min-1,  about  5  kg  of  a  catalyst  with 
0.5%  of  Pt  on  7-alumina  in  a  non-isothermal  fixed-bed  reactor  is 
needed  to  reduce  CO  from  1%  mol-basis  to  less  than  50  ppm. 

Lee  et  al.  [20]  developed  a  PrOx  reactor  for  a  lOkWe  PEMFC 
using  a  Pt-Ru/Al203 -based  catalyst.  In  order  to  diminish  the  pres¬ 
sure  drop  and  avoid  hot  spots  in  the  bed,  the  authors  proposed 
two  stage  in-series  multitubular  adiabatic  reactors  with  interme¬ 
diate  cooling.  They  reported  a  total  system  volume  about  6280  cm3, 
where  1200  cm3  correspond  to  the  catalyst  bed. 

Dudfield  et  al.  [22,23]  designed,  constructed  and  evaluated  a 
compact  unit  of  a  PrOx  reactor  treating  a  gaseous  stream  in  a 
methanol  processor  for  PEMFC  used  in  mobile  applications.  After 
evaluating  different  catalysts,  the  authors  selected  a  Pt/Ru-based 
catalyst  as  it  showed  higher  selectivity  and  conversion  on  the 
studied  temperature  range  [23].  The  system  design  consists  of  a 
Pt/Ru-based  catalyst  supported  on  the  surface  of  a  finned-type 
heat  exchanger  using  oil  as  a  coolant,  arranged  in  two  stages  in¬ 
series  of  0.5  L  each.  The  system  reduced  the  initial  CO  concentration 
from  0.5%  to  less  than  20  ppm,  operating  at  an  average  temperature 
of  160  °C,  and  considering  the  gas  flowrate  equivalent  to  produce 
5  kWe  of  electrical  power  in  a  PEMFC.  For  an  inlet  CO  composition  of 
7000  ppm  and  a  fuel  flowrate  equivalent  to  produce  3.75-7.5  kWe, 
a  unique  unit  of  0.5  L  is  not  enough  to  reduce  the  outlet  CO  level 
to  20  ppm.  Whereas  by  arranging  two  reactors  in-series  of  0.5  L 
each,  it  is  feasible  to  achieve  a  final  CO  concentration  of  10  ppm. 
By  simultaneously  increasing  the  residence  time  and  optimizing 
the  air  injection  in  each  stage,  an  exit  CO  concentration  of  10  ppm 
is  obtained  for  a  flowrate  to  produce  7.5  kWe  in  a  PEMFC.  Based  on 
the  same  concept,  Dudfield  et  al.  [22]  developed  a  reactor  to  process 
a  stream  flowrate  equivalent  to  20  kWe,  where  each  reactor  occu¬ 
pies  2000  cm3,  weighs  2.5  kg  and  has  approximately  1  L  available 
for  the  reaction  with  oxygen  and  1  L  to  maintain  the  temperature 
by  recirculation  of  a  cooling  fluid  (oil). 

Bissett  and  Oh  [  18  ]  modeled  a  monolithic-type  reactor  operating 
isothermally  and  adiabatically. 

To  the  authors’  knowledge,  model-based  optimal  reactor  design 
and  analysis  applying  mathematical  programming  techniques  have 


not  been  yet  addressed  in  this  area.  Modeling  and  optimization 
allow  evaluating  which  system  designs  and  operation  conditions 
-  e.g.  adiabatically  or  isothermally  -  are  more  appropriate  under 
different  scenarios  such  as  power  generation  scale,  costs  struc¬ 
ture,  fuel  cell  quality  requirements  (CO  content)  and  application 
purposes  (mobile  or  stationary  units). 

In  this  paper,  a  one-dimensional  heterogeneous  catalytic  reac¬ 
tor  accounting  for  interfacial  gradients  [24]  is  used  to  optimize 
the  PrOx  reactor.  Different  reactor  components  are  added  grad¬ 
ually  to  illustrate  how  the  system  dimensions  and  configuration 
changes  after  optimization.  Different  inlet  CO  compositions  and 
power  generation  targets  are  analyzed.  After  that,  simulations  with 
different  reactor  models  will  be  performed.  The  aim  of  these  sim¬ 
ulations  is  to  verify  feasibility  of  the  optimal  design  obtained  with 
the  proposed  methodology  using  the  one-dimensional  heteroge¬ 
neous  model  without  intra-particle  gradients,  which  is  intended  to 
approximate  an  egg-shell  catalyst  behavior. 


3.  Mathematical  model 

3.1.  Reactor  model 


In  this  work,  the  reactor  design  is  performed  based  on  a 
one-dimensional  heterogeneous  model.  This  model  offers  higher 
accuracy  for  reactor  design  [25].  The  mass  and  energy  balances, 
and  pressure  drop  along  the  catalytic  bed  for  the  fluid  phase  and 
balances  for  catalyst  particles  are  represented  by  the  following 
equations  [24,26]: 

Fluid  phase 

dF 

-^=Mg,iav(C,-CsSj)  (5) 

GCpr^jj  =  fifav(Tss  -  Tf)  -  Qex  (6) 


dP 

dz 


c2  1— 6b 

PfDP  eb 


+  150 


1  -  £5 

Re 


(7) 


For  solid 


Nrx 

kg.MCi  -  Cf ,.)  +  Pb^ajr/Cf,  Tf)  =  0  (8) 

j=  1 

Nrx 

hfav(Tf  -  Tl)  +  pbYJ(-AHWCl  T|)  =  0  (9) 

j=i 

Eqs.  (5)-(7)  represent  the  mass,  energy  and  momentum  balances 
for  the  fluid  phase  with  initial  conditions  at  z  =  0  given  by  F*  =  FP, 
Tf  =  Tf°  and  P=P°.  Pressure  drop  is  evaluated  using  an  Ergun-type 
expression  with  friction  parameters  given  by  Tallmadge  correla¬ 
tions  [26],  which  are  valid  for  a  wide  range  of  Reynolds  numbers 
(0.1  <Re/(  1  -£b)<  10,000).  Eqs.  (8)  and  (9)  represent  the  mass  and 
energy  balances  for  the  fluid-particle  interphase  as  boundary  con¬ 
ditions. 

3.2.  Physical- chemical  parameters  and  fluid  properties 

Ideal  gas  behavior  is  assumed  to  evaluate  the  mixture  density 
of  the  gas  phase  (pf).  The  average  molecular  weight  of  the  gaseous 
mixture  and  the  molar  flow  entering  the  reactor  were  used  to  com¬ 
pute  the  mass  flow  per  unitary  area  (G).  Correlations  for  the  fluid 
heat  capacity  (Cpr),  viscosity  (/xf)  and  thermal  conductivity  (Af)  are 
those  given  in  [27]. 
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3.2 A.  Viscosity 

The  Bromley  and  Wilke  modification  of  the  theoretical 
Hirschfleder  method  was  used  to  estimate  the  pure  component 
viscosity  [28]: 


Mi  = 


33.3 

Vc2/3 


0.261 


(1  ^(Tf/TCj))0'910810^1'9^^01^ 


(10) 


3.2.5.  Bed  properties 

The  bed  porosity  ( sb )  is  estimated  by  the  Haughey  and  Beveridge 
expression  given  by  Froment  and  Bischoff  [24]: 


£b  =  0.38  +  0.073 


(Dt/Dp  -  2 )2 
(Dt/Dp)2 


(18) 


For  spherical  shape  particles,  the  external  catalyst  surface  area 
per  unitary  reactor  volume  (av)  is  computed  by  av  =  6(l  -£b)/Dp. 
The  mass  of  catalyst  per  unitary  bed  volume  ( pb )  is  calculated 
according  to  pb  =  pp(l  -£b)* 


The  gas  mixture  viscosity  was  calculated  from  WilkeTs  method 
[27]. 


ft  f  =  1.00E-06 


y;M; 


ym 


(ii) 


where  cp^  =  <p.f  1  =  (MjAl  1)1^2  as  given  by  Fleming  and  Zipperer 
(see  Reid  et  al.  [27]) 


3.2.2.  Thermal  conductivity 

The  pure  component  thermal  conductivities  were  calculated 
using  the  Eucken’s  approximation  [28]: 

Xj  =  /U, ,4.1890  X  1(T4  (cPi 2.394  x  10"4  +  ^)  (12) 

The  following  expression  was  used  for  computing  the  mixture 
thermal  conductivity  [28,29]: 


Xf  =  0.01 


1/3 


1/3 


(13) 


3.2.3.  Heat  capacity 

The  fluid  mass  heat  capacity  is  calculated  as  cp/  = 
where  the  pure  component  heat  capacities 
(CPi)  were  obtained  from  Prausnitz  and  co-workers  [27]. 


3.2.4.  Mass  and  energy  transfer  coefficients 

The  frequently  used  correlation  for  the  mass  and  energy  trans¬ 
fer  coefficients  in  the  film,  proposed  by  Petrovic  and  Thodos  for 
3<Re<  2000  (see  Satterfield  [30]),  is  here  used. 


kg,i'Hp£b 

Difi,m 


0.357Re°641 


Mf 

PfDifjm 


1/3 


Mb 

cc« 


=  0.357 


Vi 


( 

V  Mf 


-0.359 


(  CPfMf\ 

V  ^  / 


-2/3 


(14) 

(15) 


Difz  m  is  the  diffusion  coefficient  of  component  i  in  the  gas  mixture, 
and  is  calculated  according  to  the  Wilke  equation: 


i  =  yi 

Difj  m  ^-^Di  k 

k^i 


(16) 


The  binary  diffusion  coefficients  Dik  are  given  as  a  function  of 
the  diffusion  volumes  by  the  Fuller  et  al.  estimation  (see  Reid  et  al. 
[27]): 

0.0014489475  Tf175(l/M!  +  1/Mk)1/2 

D<’k  =  ~  773  77T1  07) 

PV 2[(27v)i1/3+(2;v)y3] 


3.3.  Insulation  model 


In  Eq.  (6),  Qex  is  the  heat  exchanged  between  the  reaction  system 
and  its  surroundings.  Qex  represents  the  heat  flow  per  unitary  reac¬ 
tor  bed  volume,  and  is  related  to  the  fluid  temperature  and  the  tube 
wall  temperature  by  the  heat  transfer  coefficient  (hw)  computed  by 
the  Leva  correlation  [24]. 

Qex  =  ^w(rf-Tw)  (19) 


In  the  adiabatic  operation  mode  Qex  =  0.  In  the  thermally 
insulated  reactor  case,  a  heat  loss  towards  the  surroundings  is 
considered  by  modeling  the  energy  transfer  by  heat  conduction 
through  the  tube  and  insulation  materials,  and  by  convection  and 
radiation  from  the  insulation  surface  to  the  environment.  By  apply¬ 
ing  the  Fourier  law  to  a  cylindrical  geometry,  and  considering 
temperature  dependence  of  the  thermal  conductivity,  the  heat  flow 
by  conduction  through  the  insulation  material  results  in  the  follow¬ 
ing  expression: 


Qex 


8 

D2/n(Dins/Dt) 


k(T)dr 


(20) 


Finally,  the  heat  flow  is  related  to  the  heat  loss  from  insula¬ 
tion  surface  [31  ]  by  convection  and  radiation  by  Eqs.  (22)  and  (23), 
respectively. 


Qex  =  ^f(qcV+qrad) 


(21) 


qcv  =  2.3613  x  10 

x  (^ins  —  Tamb) 


Gcf( 


1 


ins 
1.266 


(Jamb  +  TinS)°-9  -459.67 


0.181 


(22) 


Qrad=^ins((M)4-(qmb)4)  (23) 

3.4.  Kinetics  aspects 

In  spite  of  the  large  number  of  works  proposing  kinetic  expres¬ 
sions  for  the  CO  oxidation,  only  a  few  consider  simultaneously  the 
H2  oxidation.  Other  works  report  expressions  for  the  H2  oxidation 
over  platinum  separately  [12,15-18].  Choi  and  Stenger  [5]  derived 
kinetic  expressions  for  the  simultaneous  oxidation  of  CO  and  H2  and 
the  water-gas-shift  reaction  obtained  over  a  Pt-Fe/Al203 -based 
catalyst.  As  small  numerical  values  are  computed  for  the  partial 
pressure  of  the  reactive  compounds  (02  and  CO)  it  is  necessary  to 
slightly  modify  the  original  expressions  in  order  to  avoid  numer¬ 
ical  problems  and  to  facilitate  convergence.  By  incorporating  the 
parameter  /c  =  1  x  10-6  into  Eqs.  (24)  and  (25),  the  numerical  prob¬ 
lems  are  overcome  without  loosing  accuracy  on  output  results 
(Fig.  1).  In  addition,  parameter  /c  corrects  the  reaction  kinetics  “to 
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force”  the  reaction  rates  tending  to  zero  when  the  reactant  concen¬ 
trations  tend  to  zero.  The  kinetic  expressions  are 

rco  =  0.098000  exp  ((Po2  +  k)°'5  - k03) 

(24) 

rH2  =0.005703  exp  (-^^)  ((Pq2  +  0°'5  -  fe0'5)  (25) 

rWcs  =  1.222778  exp  (-^^)  (pcoPh2o  ~  (26) 


where 


f^eq  — 


PH2PCQ2 

PcqPh2o 


being 


Fig.  2.  Catalyst  particle  models  used,  (a)  Catalyst  particle  model  used  for  optimiza¬ 
tion;  (b)  egg-shell  type  catalyst  particle  model. 

considered  above.  This  assumption  is  later  verified  in  Section 
4. 


ln(Keq)  =  56^3,5  +  1.077  x  log(T)  +  5.44  x  10-4  T  -  1.125 


x 10-7  T 2 


49170 


13.148 


3.5.  Catalyst  specifications 

The  catalyst  particle  has  a  density  of  1.1  gcat  cnrr3particle,  with  a 
Pt  density  of  0.0011  gPt  cnrr3activezone  [32],  platinum  is  the  active 
component  and  is  considered  uniformly  distributed  inside  the 
particle;  that  is,  the  whole  particle  is  assumed  as  “active  zone”. 
On  the  other  hand,  an  egg-shell  type  catalyst  particle  with  the 
same  particle  density  but  with  an  active  catalyst  density  vary¬ 
ing  along  the  radio  is  considered  (Fig.  2).  The  Pt  density  is  equal 
to  0.00771  gptcm-3activezone  in  the  “active  zone”  and  zero  in  the 
“non-active  zone”.  As  an  egg-shell  type  model  with  diffusion 
resistance  becomes  difficult  to  be  optimized  from  a  numerical 
point  of  view,  the  optimization  model  is  based  on  a  catalyst  par¬ 
ticle  with  uniform  Pt  distribution,  without  diffusion  resistance. 
This  model  is  an  approximation  to  the  egg-shell  particle  model 


3.6.  Optimization  model 


The  optimization  problem  is  formulated  to  obtain  the  optimal 
operating  conditions  and  equipment  size  aiming  at  minimizing  the 
total  system  volume.  Following,  the  optimization  problem  is  formu¬ 
lated;  i.e.  the  objective  function,  decision  variables  and  constraints 
are  specified. 


Objective  function 
Decision  variables 

Interior  point  constraints 
Catalyst  temperature 
External  insulator  temperature 

Final  point  constraints 
Admissible  pressure  drop 
CO  molar  fraction 


Min(Vtotai) 

T° ,  C,  Dt,  Dp,  eins 

j'Low  <-  ■j’S  <-  j'Up 

J.  <  jUP 
ms  — 


(Pin-Pout)/Pin<0.3 
yco  <  1  x  10-5  (10  ppm) 


Design  constraints 

Plug  flow  condition  Lt /Dp  >  50;  Dt/Dp  >10 

As  mentioned,  it  is  intended  to  minimize  the  total  volume  of 
the  system.  In  a  conventional  fixed-bed  reactor,  the  total  volume  is 
computed  considering  the  different  reactor  items  such  as  the  cat¬ 
alyst  bed,  reactor  tube,  insulating  material,  distribution  heads  and 
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entrance  length  with  inert  material,  which  determine  the  following 
individual  volumes: 

reactor  volume  :  Vreactor  =  ^t  +  2^ns) 

.  .  ,  ,  20Dp(Dt  +  2ems)2n 

inert  material  volume  :  Vinert  material  = - - - ^ - 

hemispherical  heads  volume  : 

x/  _  4(Dt  +  2eins/2)37 r 

v hemispherical  heads  —  ^ 

Then,  the  total  volume  is 


Hemispherical  Heads 
Distributor  -  Collector 


^total  —  ^reactor  +  ^hemispherical  heads  +  ^inert  material 

The  optimization  problem  determines  the  optimal  reactor 
length  (Lt),  reactor  diameter  (Dt),  catalyst  particle  diameter  (Dp) 
and  insulating  material  thickness  (eins)  that  minimize  the  total 
system  volume.  For  adiabatic  operation  case,  the  constraints  cor¬ 
responding  to  the  insulating  material  temperature  (Tins)  and 
insulation  thickness  are  not  included  into  the  model  as  decision 
variables. 

Due  to  the  exothermic  nature  of  the  PrOx  reaction,  it  is  desir¬ 
able  to  keep  the  reactor  temperature  within  a  certain  operating 
range.  An  upper  temperature  bound  is  set  for  avoiding  catalyst  sin¬ 
tering.  At  low  temperature,  the  reaction  rate  diminishes,  the  kinetic 
expression  is  not  longer  valid,  and  water  condensation  occurs.  Then, 
it  is  necessary  to  impose  a  lower  temperature  bound.  The  reactor 
inlet  temperature  (T°)  is  a  decision  variable,  and  its  value  results 
from  the  optimization  problem. 

Summarizing,  molar  flow  of  component  i  (F{),  fluid  tempera¬ 
ture  (Tf)  and  pressure  along  the  reactor  bed  (P)  are  calculated  from 
differential  Eqs.  (5)-(7).  Component  concentration  (C|  f)  and  tem¬ 
perature  (Tf)  at  the  external  catalyst  surface  are  obtained  from  Eqs. 
(8)  and  (9).  The  resulting  differential  algebraic  equations  (DAEs)  are 
implemented  and  solved  using  gPROMS  (general  process  modeling 
system)  [33].  gPROMS  is  a  general  purpose  modeling,  simulation 
and  optimization  system  software.  DAEs  are  essentially  sets  of  ordi¬ 
nary  differential  equations  (ODEs)  where  some  of  the  variables  are 
constrained  by  algebraic  relations.  Then,  the  resulting  initial  value 
DAE  problem  is  integrated  along  the  z  axial  direction.  Backward 
differentiation  formula  (BDF)  type  method  is  the  algorithm  chosen 
in  gPROMS  for  solving  the  differential  equation  system.  Finally,  the 
optimization  algorithm  being  used  is  the  single-shooting  method, 
which  is  also  available  in  the  gPROMS  environment. 

4.  Results  and  discussion 

The  reactor  is  designed  considering  the  most  probable  condi¬ 
tions  for  the  PrOx  reaction  in  a  small  or  medium  scale  ethanol 
processor.  Molar  flow  rates  and  compositions  at  the  PrOx  reactor 
inlet  considered  for  design  correspond  to  the  WGS  exit  conditions 


analyzed  in  a  previous  work  [34].  The  inlet  molar  flow  rates  per  kW 
considered  at  the  design  stage  are  shown  in  Table  1 .  Apart  from  the 
main  process  stream,  oxygen  is  fed  to  the  reactor  to  perform  the 
oxidation  reaction.  The  results  presented  in  this  work  correspond 
to  an  inlet  02-to-CO  molar  flowrate  ratio  equal  to  2. 

4.1.  Design  for  0.3%  CO  at  the  reactor  inlet 

Firstly,  a  system  consisting  on  a  single-stage  reactor  is  con¬ 
sidered.  Four  different  objective  functions  are  formulated,  which 
differ  in  the  number  of  system  components  included  into  the 
unit  size  computation.  Specifically,  Table  2  shows  the  results 
obtained  for:  (a)  adiabatic  reactor,  (b)  reactor  +  insulation,  (c)  reac¬ 
tor  +  insulation  +  distribution  heads  and  finally,  (d)  same  as  case  (c) 
but  evaluating  also  an  inert  entrance  length.  A  0.3%  CO  composition 
at  the  reactor  inlet  is  assumed  for  all  cases. 

In  the  adiabatic  case,  only  the  volume  of  the  catalyst  bed  is  con¬ 
sidered.  The  “reactor  +  insulation  material”  system  considers  the 
reactor  tube  and  insulation  volumes.  The  insulating  material  used 
in  the  model  consists  of  a  layer  of  calcium  silicate  covered  with 
a  thin  layer  of  aluminum.  The  insulating  material  conductivity  is 
considered  as  a  quadratic  function  of  temperature  according  to  the 
correlation  k(T)  =  0.0264  +  7  x  10_5T  +  6  x  10_8T2  taken  from  Indus¬ 
trial  Insulation  Group  brochure  [35]. 

A  fully  developed  flow  is  considered  at  the  reactor  bed  entrance. 
However,  it  should  be  taken  into  account  that  the  flow  into  the 
fixed-bed  reactor  is  generally  achieved  by  means  of  a  feed  pipe  and 
a  distribution  hood.  These  must  therefore  be  constructed  so  that 
the  fixed  bed  is  uniformly  traversed,  and  the  gas  residence  time  in 
each  flow  path  of  the  fixed  bed  is  the  same.  In  order  to  evaluate  the 
relative  size  of  the  components  to  achieve  a  uniform  gas  flow  to  the 
catalyst  bed,  the  model  is  refined  by  including  the  reactor  heads 
and  an  entrance  with  spherical  inert  particles.  It  has  been  found 
that  an  entrance  length  about  Le  =  20DP  is  enough  to  approach  con¬ 
stant  values  independent  of  I  for  the  gas  flow  parameters  [29]. 
Finally,  hemispherical  distributor  and  collector  are  considered  at 
the  beginning  and  ending  of  the  reactor  (Fig.  3). 

Table  2  lists  the  optimization  results  for  the  four  cases  analyzed 
to  produce  lOkW  for  a  0.3%  CO  inlet  concentration.  Evidently,  the 


Table  1 

Inlet  PrOx  reactor  molar  flow  rates  and  compositions 


Design 

1%C0 

0.7%  CO 

0.3%  CO 

Molar  flow  rate  (mol  h_1 ) 

Molar  fraction 

Molar  flow  rate  (mol  h-1 ) 

Molar  fraction 

Molar  flow  rate  (mol  hr1 ) 

Molar  fraction 

ch4 

4.67 

0.0563 

4.67 

0.0563 

4.67 

0.0563 

h2 

41.97 

0.5057 

42.22 

0.5087 

42.55 

0.5127 

CO 

0.83 

0.0100 

0.58 

0.0070 

0.25 

0.0030 

h2o 

20.59 

0.2481 

20.34 

0.2451 

20.01 

0.2411 

co2 

14.94 

0.1800 

15.19 

0.1830 

15.52 

0.1870 
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Table  2 

Optimization  results  considering  different  objective  functions 
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Design  for  0.3%  CO 

Power  target  10  kW 

(a) 

(b) 

(c) 

(d) 

U  (cm) 

2.81 

2.5* 

9.37 

8.25 

Dt  (cm) 

6.19 

6.74 

3.71 

4.08 

Dp  (cm) 

0.05* 

0.05* 

0.19 

0.10 

eais  (cm) 

0 

3.51 

2.91 

2.91 

Tin  (K) 

528.0 

522.0 

511.4 

507.6 

Tout  (K) 

574.0 

567.6 

554.8 

550.2 

Bed  volume  (cm3) 

84.5(100%) 

89(24%) 

101 (09%) 

108(08%) 

Insulation  volume  (cm3) 

- 

282(76%) 

567(51%) 

527(41%) 

Heads  volume  (cm3) 

- 

- 

452(40%) 

508(39%) 

Inert  entrance  volume  (cm3) 

- 

- 

- 

154(12%) 

Total  volume  (cm3) 

84.5 

372 

1120 

1297 

Bed  volume  per  kW  produced  (cm3  1<W-1 ) 

8.5 

8.9 

10.13 

10.79 

Insulation  volume  per  kW  produced  (cm3  1<W-1 ) 

0 

28.2 

56.66 

52.72 

Heads  volume  per  kW  produced  (cm3  kW-1 ) 

0 

0 

45.24 

50.8 

Entrance  volume  per  kW  produced  (cm3  1<W-1 ) 

0 

0 

0 

15.4 

Total  volume  per  kW  produced  (cm3  1<W-1 ) 

8.5 

37.2 

112.0 

129.7 

(a)  Adiabatic  reactor;  (b)  reactor  +  insulation;  (c)  reactor  +  insulation  +  heads;  (d)  reactor  +  insulation  +  heads  +  inert  entrance;  the  asterisk  shows  that  the  variable  reaches  the 
minimum  admitted  value  in  the  optimization  process. 


more  detailed  the  model  is,  the  better  it  represents  the  real  sys¬ 
tem  investigated.  The  adiabatic  reactor  model  computes  an  optimal 
reactor  diameter  greater  than  its  length,  rendering  a  diameter-to- 
length  ratio  of  2.2.  Then,  the  model  is  extended  by  adding  insulation 
material.  In  this  case,  the  catalytic  bed  volume  increases  5.3%  with 
respect  to  the  previous  model,  and  the  diameter-to-length  ratio 
increases  to  2.7.  It  should  be  noted  that  the  computed  reactor  length 
is  at  the  lower  bound  value  fixed.  This  is  because  the  model  tries 
to  take  advantage  of  the  exothermic  reaction  which  is  developed 
inside  the  reactor  along  radial  direction.  The  total  system  volume 
is  372  cm3,  of  which  76%  corresponds  to  the  insulating  material 
volume.  If  hemispherical  heads  -  to  favor  an  uniformly  distributed 
flux  -  are  incorporated  into  the  model,  the  optimization  results 
provide  a  different  “qualitative”  reactor  geometry;  specifically,  in 
this  case  the  optimal  reactor  length  is  greater  than  its  diameter 
instead,  rendering  a  diameter-to-length  ratio  of  0.4.  The  total  sys¬ 
tem  volume  (1120  cm3 )  increases  tremendously  with  respect  to  the 
previous  cases.  The  volume  of  the  catalytic  bed,  insulating  material 
and  heads  represent  9%,  51%  and  40%  of  the  total  volume,  respec¬ 
tively.  Finally,  an  inert  length  entrance  is  also  incorporated  into  the 
model  for  the  same  reason  as  reactor  heads  were  considered.  In 
this  case,  the  diameter-to-length  ratio  is  around  0.5.  The  optimal 
volume  computed  for  the  entrance  length  is  154  cm3,  which  rep¬ 
resents  12%  of  the  total  system  volume.  The  total  volume  increases 
16%  with  respect  to  the  previous  case.  Based  on  these  model-based 
results,  the  final  reactor  design  is  mainly  governed  by  the  presence 
of  hemispherical  heads. 

Table  2  shows  that  the  inlet  reactor  temperature  Tin  varies  20  K 
from  model  (a)  -  column  1  -  to  model  (c)  -  column  4  -  and  that 
the  temperature  differences  at  the  reactor  outlet  (Tou t)  and  inlet 
(Tin)  range  between  42.6  K  and  46  K  for  all  analyzed  cases.  The 
small  disagreement  among  models  is  because  of  the  variation  of 
the  reaction  heat  with  temperature  and  the  heat  loss  to  surround¬ 
ings. 

Finally,  it  can  be  observed  that  the  optimal  catalyst  particle  size 
is  sensitive  to  the  model  being  used.  Certainly,  the  particle  diameter 
computed  for  models  (a)  and  (b)  is  at  the  lower  bound  value  fixed 
(0.05  cm),  and  0.19  and  0.10  cm  for  models  (c)  and  (d),  respectively. 
A  lower  bound  for  particle  diameter  of  0.05  cm  was  set  for  practical 
reasons.  In  order  to  compare  optimization  results  with  respect  to 
other  scenarios,  Table  2  also  lists  the  optimal  volume  values  per  kW 
of  generated  power. 


4.2.  Effect  of  the  power  generation  level 

Four  design  targets  regarding  the  fuel  cell  power  generation 
level  are  analyzed  to  achieve  10  ppm  CO  concentration  at  the  reac¬ 
tor  outlet:  1, 10,  25  and  50  kW.  As  above,  process  streams  with  CO 
compositions  of  0.3%,  0.7%  and  1%  at  the  reactor  inlet  are  assumed. 
The  results  obtained  are  listed  in  Table  3.  The  optimization  prob¬ 
lem  allowed  achieving  a  feasible  design  for  each  case  that  fulfills 
the  imposed  constraints.  It  should  be  noted  that  different  system 
designs  and  arrangements  are  achieved  depending  on  the  desired 
power  targets. 

As  mentioned,  due  to  the  exothermic  nature  of  the  involved 
reactions,  a  main  difficulty  in  reactor  design  is  to  keep  the  operation 
temperature  within  the  desired  limits.  If  1%  CO  level  is  specified 
at  the  reactor  inlet,  two  units  in-series  (R1  and  R2)  are  needed 
to  achieve  a  CO  concentration  of  10  ppm  at  the  reactor  outlet  for 
the  four  power  generation  targets  being  considered.  That  is,  an 
arrangement  of  two  units  with  intermediate  cooling  has  to  be 
considered  for  those  scenarios.  Figs.  4-6  show  the  relative  vol¬ 
umes  (in  percentage)  computed  for  the  system  components  for 
all  power  design  targets  and  inlet  CO  compositions  being  inves¬ 
tigated.  Since  the  reaction  is  highly  exothermic,  a  large  insulation 
volume  is  required.  Heads  and  entrance  lengths  have  an  important 
contribution  to  the  total  volume  of  the  system.  That  is  necessary 
to  achieve  a  uniform  gas  flow  distribution  pattern  to  the  catalyst 
bed. 


kW  generated 


Fig.  4.  Relative  volumes  of  system  components  corresponding  to  four  generation 
power  targets  for  0.3%  CO  molar  fraction  at  the  entrance. 
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Table  3 

Optimization  results  for  different  design  targets 


Reactor  +  Insulation  +  Headers  +  Entrance  Length 

CO  Inlet  [%] 

0.3 

0.7 

1 

CO  Outlet  [ppm] 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

Power  [kW] 

1 

10 

25 

50 

1 

10 

25 

50 

1 

10 

25 

50 

Amount  of  units 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

Unit# 

R1 

R1 

R1 

R1 

R1 

R1 

R1 

R1 

R1 

R2 

R1 

R2 

R1 

R2 

R1 

R2 

Lt  [cm] 

6.08 

8.25 

8.45 

8.94 

6.50 

9.88 

11.54 

13.49 

4.77 

5.50 

9.53 

5.84 

10.51 

4.15 

6.20 

6.20 

Dt  [cm] 

1.85 

4.08 

6.21 

8.20 

1.88 

4.89 

7.19 

9.45 

2.08 

1.58 

4.78 

4.53 

8.50 

7.08 

9.72 

9.52 

Dp  [cm] 

0.05 

0.10 

0.12 

0.16 

0.05 

0.06 

0.09 

0.14 

0.05 

0.05 

0.12 

0.09 

0.07 

0.08 

0.12 

0.12 

eins  [cm] 

1.97 

2.91 

3.32 

3.69 

2  56 

3.36 

3.71 

3.98 

2.06 

2.43 

2.76 

3.24 

3.85 

2.34 

3.97 

4.00 

Tin  [K] 

474.4 

507.6 

511.8 

518.7 

475.3 

475.3 

474.0 

472.8 

451.4 

499.0 

452.5 

505.0 

447.7 

465.7 

497.9 

508.7 

T  out  [K] 

516.8 

550.2 

555.3 

562.4 

571.8 

573.9 

573.5 

572.7 

519.9 

568.6 

531.9 

567.8 

574.3 

485.7 

570.2 

572.2 

Bed  Vol.  [cmJ] 

16 

108 

256 

472 

18 

186 

469 

946 

16 

11 

171 

94 

596 

163 

460 

441 

Insul.Vol.  [cm'3] 

144 

527 

840 

1232 

232 

860 

1466 

2265 

128 

168 

623 

462 

1570 

287 

1059 

1053 

Heads  Vol.  [cm3] 

102 

508 

1111 

1980 

180 

819 

1633 

2763 

125 

140 

572 

699 

2226 

852 

2884 

2816 

Entrance  Vol.  [cm3] 

26 

154 

311 

610 

38 

127 

302 

667 

30 

33 

200 

171 

289 

174 

588 

579 

Total  Vol.  [cm3] 

288 

1297 

2518 

4295 

468 

1992 

3869 

6641 

651 

2992 

6157 

9879 

Bed  Vol./kW  [cmJ 
kW1] 

16.35 

10.79 

10.24 

9.44 

18.04 

18.56 

18.74 

18.92 

16.21 

10.78 

17.10 

9.41 

23.86 

6.54 

9.20 

8.83 

Insul  Vol./kW  [cm* 
kW’l 

143.84 

52.72 

33.60 

24.64 

232.11 

86.04 

58.64 

45.31 

127.80 

168.37 

62.31 

46.19 

62.80 

11.50 

21.17 

21.07 

Heads  Vol./kW  [cm3 
kW’l 

101.63 

50.80 

44.44 

39.60 

179.59 

81.94 

65.31 

55.26 

124.79 

139.85 

57.22 

69.88 

89.04 

34.06 

57.68 

56.32 

Entrance  Vol./kW 
[cm3  kW1] 

26.33 

15.40 

12.45 

12.20 

38.48 

12.70 

12.07 

13.33 

30.19 

32.57 

20.00 

17.14 

11.54 

6.95 

11.76 

11.57 

Total  Vol./kW  [cm3 
kW1! 

288 

130 

101 

86 

468 

199 

155 

133 

651 

299 
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198 
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Fig.  5.  Relative  volumes  of  system  components  corresponding  to  four  generation 
power  targets  for  0.7%  CO  molar  fraction  at  the  entrance. 
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Fig.  6.  Relative  volumes  of  system  components  corresponding  to  four  generation 
power  targets  for  1%  CO  molar  fraction  at  the  entrance.  Two  reaction  units  in-series 
(R1  and  R2). 


4.3.  Reactor  simulation 

Following,  a  comparison  of  simulated  profiles  obtained  with 
different  reactor  models  is  presented  and  discussed.  In  order  to 
verify  feasibility  of  the  optimal  design  obtained  with  the  method¬ 
ology  proposed  using  the  one-dimensional  heterogeneous  model 
without  intra-particle  gradients,  which  is  intended  to  approximate 
the  egg-shell  catalyst  behavior,  simulation  runs  were  then  con¬ 
ducted  based  on  different  reactor  models.  The  simulations  results 
correspond  to  the  scenario  with  an  inlet  CO  concentration  of 
0.7%  and  a  flow  rate  equivalent  to  generate  10  kW  of  power  (see 
Table  3). 

The  reactor  models  being  compared  are:  (i)  the  mentioned 
model  used  for  optimal  design,  named  lD-HtEx;  (ii)  a  one- 
dimension  heterogeneous  model  with  intra-particle  gradients, 
named  lD-Htlnt;  finally,  (iii)  a  two-dimensional  pseudohomoge- 
nous  model,  named  2D-PH.  Based  on  the  lD-Htlnt  model, 
two  different  situations  regarding  the  catalyst  distribution  were 
addressed:  (a)  uniform  distribution  of  the  catalyst  activity  in  the 
particle  with  a  Pt  density  equal  to  0.0011  gptcnrr  3  cat.  named  1D- 
HtlntU;  (b)  an  “egg-shell”-type  model,  named  lD-HtlntES,  where 
the  catalyst  activity  is  concentrated  in  an  “active  outer  zone”  corre¬ 
sponding  to  around  5%  of  the  particle  radius  with  a  Pt  density  equal 
tO  0.00771  gpt  Cm  ^active  zone  [32]. 

The  CO  molar  fraction  profiles  established  along  the  reactor  are 
depicted  in  Fig.  7,  and  the  composition  profiles  inside  the  catalyst 
particle  are  shown  in  Fig.  8.  It  can  be  observed  that  the  conversion 
values  predicted  by  the  lD-HtEx  model  and  the  lD-HtlntES  model 
are  similar.  Meanwhile  the  2D-PEI  model  and  1  D-EItlntU  model  pre¬ 
dict  higher  and  lower  conversion  values  than  the  lD-EItEx  model, 
respectively. 

This  behavior  is  due  to  the  fact  that  the  average  temperature 
predicted  by  the  2D-PEI  model  is  higher  (Fig.  9),  which  increases 
the  reaction  rate,  thus  requiring  less  volume  to  achieve  the  desired 
target.  Although  the  conversion  predicted  by  the  lD-EItlntU  model 
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Dimensionless  catalyst  particle  radius 


Fig.  8.  CO  molar  fraction  profiles  inside  the  catalyst  particle. 


Dimensionless  axial  coordinate 

Fig.  9.  Temperature  profiles  predicted  by  the  different  models  considered. 


is  higher  at  the  first  sections  of  the  reactor,  which  is  indicated 
through  the  higher  temperature  profile,  the  diffusion  phenomena 
compete  with  the  reaction  ones  at  the  last  sections  of  the  unit,  thus 
decreasing  the  catalyst  effectiveness  and,  consequently,  lowering 
the  conversion  values. 

According  to  the  aims  set  out  in  this  paper,  the  results  of  the 
models  comparison  allow  confirming  the  application  of  the  1D- 
HtEx  model  for  optimal  reactor  design  when  considering  a  catalyst 
with  non-uniformly  distributed  activity  (egg-shell  type  catalyst). 
Nevertheless,  the  higher  average  temperature  predicted  by  the  2D- 
PH  model  suggests  further  consideration  of  a  two-dimensional 
heterogeneous  model.  In  addition,  in  future  works,  the  02-to-CO 
molar  ratio  will  be  considered  as  a  decision  variable,  as  well  as 
the  possibility  of  oxygen  and/or  water  flow  distribution  along  the 
reaction  zone  instead  of  a  unique  feed  point. 


5.  Conclusions 

As  regards  the  used  methodology,  the  results  of  the  present 
work  reflect  clearly  the  advantages  of  applying  mathematical 
programming  techniques  to  optimize  both  design  and  operation 
conditions  of  the  PrOx  reactor,  which  is  subjected  to  several  trade¬ 
offs  involving  operative,  construction,  technological  and  efficiency 
constraints.  Indeed,  the  proposed  methodology  was  applied  to  dif¬ 
ferent  scenarios  and  targets,  showing  satisfactory  robustness  and 
flexibility. 

The  heterogeneous  reactor  model  being  used  allows  computing 
the  optimal  reactor  length  and  diameter,  optimal  catalyst  parti¬ 
cle  diameter,  optimal  insulating  material  thickness,  as  well  as  the 
optimal  inlet  temperature  of  the  stream  to  operate  the  system  in  a 
pseudo-adiabatic  mode. 

Under  the  investigated  scenarios,  and  depending  on  the  CO  con¬ 
tent  to  be  treated,  more  than  one  catalytic  bed  may  be  required. 
Indeed,  it  is  possible  to  successfully  process  a  stream  with  an  inlet 
CO  concentration  of  0.3%  or  0.7%  in  a  unique  stage,  maintaining 
restrictions  on  the  bed  temperature.  However,  for  higher  CO  con¬ 
centration,  1%  for  instance,  two  stages  in-series  with  intermediate 
cooling  are  needed. 

The  model  being  used  predicts  a  behavior  similar  to  a  catalyst 
with  a  non-uniformly  distributed  activity  (“egg-shell”).  Neverthe¬ 
less,  the  2D-PH  model  and  lD-HtlntU  model  predict  higher  and 
lower  conversions  with  respect  to  the  ones  corresponding  to  the 
design  targets. 

Although  the  presented  optimization  results  strongly  depend 
on  the  particular  catalyst  and  on  the  considered  inlet/outlet  speci¬ 
fications,  the  proposed  methodology  and  the  obtained  preliminary 
results  can  assist  in  designing,  optimizing  and  controlling  the  global 
process,  i.e.  commercial  fuel  processors  for  producing  hydrogen  for 
fuel  cells  applications. 
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